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Abstract. We investigate the possibility of the r-process during the magnetohydrohynamical ex- 
plosion of supernova in a massive star of 13 solar mass with the effects of neutrinos induced. We 
adopt five kinds of initial models which include properties of rotation and the toroidal component 
of the magnetic field . The simulations which succeed the explosions are limitted to a concentrated 
magnetic field and strong differential rotation. Low Y e ejecta produce heavy elements and the third 
peak can be reprocuced. However, the second peak is low because Y e distribution as a function of 
radius is steep and ejecta corresponding to middle Y e is very few. 
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INTRODUCTION 

It has been considered that the origin of heavy neutron-rich elements like uranium is 
mainly due to the r-process nucleosynthesis that occurs during the supernova explosions 
and/or neutron star mergers [1, 2, 3, 4]. The main issue concerning the r-process research 
is to reproduce the three peaks (A ~ 80, 130, and 195) in the abundance pattern for the 
r-elements in the solar system. Among models of the r-process, it has been believed that 
supernovae are the most plausible astrophysical site [5]. 

In the present paper, we will carry out the calculations of the MHD explosion of the 
He-core of 3.3 M during tf ~ 500 ms whose mass in the main sequence stage is about 
13 Mq star. For the MHD calculations, five distributions of initial rotation and magnetic 



TABLE 1. Initial parameters of precollapse models. 





T/\W\ (%) 


E m /\W\ (%) 


X (10 8 cm) 


Z (10 8 cm) 


n (s- 1 ) 


Bo (G) 


model 1 




















model 2 


0.5 


0.1 


1 


1 


5.2 


5.4 xlO 12 


model 3 


0.5 


0.1 


0.5 


1 


7.9 


1.0 xlO 13 


model 4 


0.5 


0.1 


0.1 


1 


42.9 


5.2 xlO 13 


model 5 


1.5 


0.1 


0.1 


1 


72.9 


5.2 xlO 13 



fields are assumed with parametric forms. We include the effects of neutrinos by using 
Leakage scheme [6]. Thereafter, we investigate the r-process nucleosynthesis using the 
results of MHD calculations. 



INITIAL MODELS 



The presupernova model has been calculated from the evolution of He-core of 3.3 
M Q that corresponds to 13 M & in the main sequence stage [7]. We adopt cylindrical 
properties of the angular velocity Q. and the toroidal component of the magnetic field 
as follows [10]: 



Q^ZH^xi^, W)=So x (1) 

where X and Z are the distances from the rotational axis and the equatorial plane with 
Xo and Zq being model parameters. 

Both Q.0 and Bo are the initial values at X = and Z = 0. Initial parameters of five 
precollapse models are given in Table 1 . The spherically symmetric case is denoted by 
model 1 . In model 2, the profiles of rotation and magnetic field in the Fe-core are taken 
to be nearly uniform. We present model 4 and 5 as the case having a differentially rapid 
rotating core and strong magnetic fields. An intermediate example, model 3 between 
model 2 and model 4 is prepared for reference. 



EXPLOSION MODELS AND r-PROCESS NUCLEOSYNTHESIS 

We perform the calculations of the collapse, bounce, and the propagation of the shock 
wave with use of ZEUS-2D [8] in which the realistic equation of state [9] has been 
implemented by Kotake et al[10]. It is noted that the contribution of the nuclear energy 
generation is usually negligible compared to the shock energy. In Table 2, our results of 
MHD calculations are summarized. 

The simulations which succeed the explosions are model 4 and model 5. Explosion 
fails for spherical (modell) and/or large Xq models, which indicate that explosions need 
concentrated magnetic field and strong differential rotation. Model 4 eject large area 
around rotation axis, but this area has high Y e . Theafore elements responsible for the 
third peak do not appear. On the other hand, Model 5 ejects collimated jet and this 
ejecta from deep area has very low Y e . Therefore we investigate the possibility of the 
r-process for model 5. Low Y e ejecta produces heavy elements and the third peak can 
be reprocuced(Fig. 2). However, the second peak is low because Y e distribution as a 
function of radius is steep and ejecta corresponding to middle Y e is very few(Fig. 1). 



SUMMARY 



We calculated magneto-hydrodynamical simulations with neutrino effects included and 
investigated the r-process nucleosynthesis during the purely magnetohydrodynamic ex- 



TABLE 2. Hydrodynamics simulation results. 





tb 


T/\W\ f 


Em/\W\f 


TP* 
^exp 


M ej /M 


M re j/M 


model 1 


111 


283 








0.023 




model 2 


125 


311 


6.91 


0.053 


0.127 




model 3 


129 


329 


8.74 


0.116 


0.164 




model 4 


133 


433 


8.80 


0.142 


1.13 


0.111 


model 5 


180 


548 


15.3 


0.339 


0.484 


0.022 



% indicates the time (ms) at the bounce. The calculations are stopped at the time tt (ms). The ratios 
T/\W\f and E m /\W\f are expressed in %. E* xp = £ e xp/10 51 ergs. 



1.0-10 



model4 
1 mode!5 



E 1.0-10" 2 



1.0-10' 



1.0-10" 



1.0-10"- 



0.05 0.1 



t 10 " 6 i % 



10 c 



solar: 
mode!5: 



10" 10 
10''" 



150 
mass number 



0.2 0.3 0.4 0.5 50 100 

Ye 

FIGURE 1. Ejected mass as a function of Y e FIGURE 2. Abundances obtained from model5. 
(model4 and model5). 



plosion in a massive star of 13 M Q progenitor model by [7]. Weak differential rotation 
model does not succeed in explosion. So, explosion need very strong and concentrated 
rotation, and magnetic field. Low Y e ejecta from which the third peak elements appear 
needs strong differential rotation and concentrated magnetic field. Explosion depend on 
the parameter of magnetic field and rotation profile and Y e distribution of ejecta is very 
different, which depend on each parameter. For model 5, third peak can be reprocuced 
but the second peak is low because Y e distribution as a function of radius is steep. 
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